Experimental vs. Numerical Eigenvalues of a Bunimovich Stadium Billiard 

A Comparison 
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We compare the statistical properties of eigenvalue sequences for a 7 = 1 Bunimovich stadium 
billiard. The eigenvalues have been obtained by two ways: one set results from a measurement of 
the eigenfrequencies of a superconducting microwave resonator (real system) and the other set is 
calculated numerically {ideal system). The influence of the mechanical imperfections of the real sys- 
tem in the analysis of the spectral fluctuations and in the length spectra compared to the exact data 
of the ideal system are shown. We also discuss the influence of a family of marginally stable orbits, 
the bouncing ball orbits, in two microwave stadium billiards with different geometrical dimensions. 

PACS number(s); 2.60.Cb, 03.65. Ge, 05.45.Mt, 41.20.Cv 



I. INTRODUCTION 



Quantum manifestations of classical chaos have re- 
ceived much attention in recent years B and for the 
semiclassical quantization of conservative chaotic sys- 
tems, two-dimensional biUiard systems provide a very 
effective tool |§,||. Due to the conserved energy of the 
ideal particle propagating inside the billiard's boundaries 
with specular reflections on the walls, the billiards be- 
long to the class of Hamiltonian systems with the low- 
est degree of freedom in which chaos can occur and this 
only depends on the given boundary shape. Such sys- 
tems are in particular adequate to study the behavior of 
the particle in the corresponding quantum regime, where 
spectral properties are completely described by the sta- 
tionary Schrodinger equation. The spectral fluctuations 
properties of such systems were investigated both ana- 
lytically and numerically. It has been found that these 
properties coincide with those of the ensembles of ran- 
dom matrix theory (RMT) having the proper symmetry 
|Q,g| if the given system is classically non-integrable. For 
time-reversal invariant systems, to which group the here 
investigated billiards belong, the relevant ensemble is the 
Gaussian orthogonal ensemble (GOE). 

In the last decades this subject was dominated by the- 
ory and numerical simulations. About seven years ago 
experimentalists have found effective techniques to sim- 
ulate quantum billiard problems with the help of macro- 
scopic devices. Due to the equivalence of the stationary 
Schrodinger equation and the classical Helmholtz equa- 
tion in two dimensions one is able to model the billiard 
by a similarly shaped electromagnetic resonator 
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FIG. 1. Typical transmission spectrum of the supercon- 
ducting 7 = 1 stadium billiard fabricated from niobium in 
the range between 10.5 and 11.5 GHz taken at 4.2 K. The 
signal is given as the ratio of output power to input power 
on a logarithmic scale. The inset illustrates the shape of the 
resonator and the positions of the antennas. 

Theoretical predictions assume to have always an ideal 
system with a perfect geometry whereas experiments 
have been performed with real systems. Real microwave 
cavities - in particular the here used superconducting 
ones - are usually not machined by the most accurate 
technique, e.g. milling from a solid block. They are cut 
from niobium sheet material which is shaped and after- 
wards electron beam welded. Finally they are chemically 
etched to clean their inner surface. Their final inner ge- 
ometry is not accessible for a direct measurement so that 
their exact geometric properties are not known. In the 
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case of the investigated Bunimovich stadium billiard [^ , 
where 7 = a/R = 1 (see inset of Fig. 1), the following 
properties are especially crucial: The radius of curvature 
of the boundary does not change abruptly but smoothly 
at the transition from the straight a to the circular sec- 
tion of the boundary |1^ . Another point is that the angle 
at the corners is not exact 90 degrees. 

The paper is about the question to which extent a com- 
parison of experimental data with theoretical predictions 
for such billiards is meaningful. Therefore we want to 
compare a numerical simulation (calculation of eigenval- 
ues) for a 7 = 1 stadium billiard with a measurement of 
a real superconducting microwave cavity (measurement 
of eigenvalues) by studying the statistical properties of 
the two sequences of eigenvalues. 

The paper is organized as follows. In Sec. || the ex- 
perimental set-up and the measurement of the eigenfre- 
quencies are described and in Sec. Ill the numerical cal- 
culations. Th e co mparison of both sets of eigenvalues is 
shown in Sec. IV by analyzing their spectral fluctuations. 



II. EXPERIMENT 

Experimentally we have investigated a two- 
dimensional microwave cavity which simulates a two- 
dimensional quantum billiard. Here we present results 
based on measurements using a superconducting niobium 
cavity, having the shape of a quarter Bunimovich stadium 
billiard. The billiard has been desymmetrized to avoid 
superpositions of several independent symmetry classes 
PI . Its inner dimensions are a = i? = 20 cm correspond- 
ing to 7 = a/R = 1, see inset of Fig. 1, and it has a height 
of d = 0.7 cm, that guarantees a two-dimensionality up 
to a frequency of fmax = c/2d = 21.4 GHz (c denotes the 
speed of light). 

As in previous investigations [n^JlSl the measurement 
of the 7=1 stadium billiard has been carried out in 
a LHe-bath cryostat. This experimental set-up is very 
stable concerning temperature and pressure fluctuations. 
The cavity has been put into a copper box which was 
covered by liquid helium, so that a constant temperature 
of 4.2 K inside the resonator was guaranteed during the 
whole measurement. The box has also been evacuated to 
a pressure of 10~^ mbar to eliminate effects of the dielec- 
tric gas inside the cavity. We were able to excite the cav- 
ity in the frequency range of 45 MHz < / < 20 GHz in 10 
kHz-steps using four capacitively coupling dipole anten- 
nas sitting in small holes on the niobium surface, see inset 
of Fig. 1 . These antennas penetrated only up to a maxi- 
mum of 0.5 mm into the cavity to avoid perturbations of 
the electromagnetic field inside the resonator. Using one 
antenna for the excitation and either another one or the 
same one for the detection of the microwave signal, we are 
able to measure the transmission as well as the reflection 
spectra of the resonator using an HP8510B vector net- 
work analyzer. In Fig. 1 a typical transmission spectrum 



of the billiard in the range between 10.5 and 11.5 GHz is 
shown. The signal is given as the ratio of output power to 
input power on a logarithmic scale. The measured reso- 
nances have quality factors of up to Q = //A/ « 10^ and 
signal-to-noise ratios of up to S/N k, 60 dB which made 
it easy to separate the resonances from each other and 
detect also weak ones above the background. As a conse- 
quence of using superconducting resonators, all the im- 
portant characteristics like eigenfrequencies and widths 
can be extracted with a very high accuracy |lj-|l^] . A de- 
tailed analysis of the raw spectra yielded a total number 
of 955 resonances up to 20 GHz. To reduce the possibil- 
ity of missing certain modes with a rather weak electric 
held vector at the position of the antennas, the measure- 
ments were always performed with different combinations 
of antennas. Thereby the number of missed modes is 
dramatically reduced below three to flve in a typical case 
of a measurement of about one thousand eigenfrequen- 
cies. For the 7=1 stadium billiard investigated here 
the 955 detected eigenmodes agrees exactly with the ex- 
pected number calculated from particular geometry of 
the cavity with the help of Weyl's formula (see Sec. IV 
below) . The measured sequence of frequencies extracted 
from the experimental spectra forms the basic set of the 
statistical investigations in Sec. IV. 



III. NUMERICAL CALCULATIONS 

A. Theory 

The problem that we have to solve is to flnd the eigen- 
values of the Dirichlet problem in a billiard, i.e. to solve 
the following equations: 



A*(f) + k^^{r) = for f e P 
^(f) =0 for f Gap 



(1) 
(2) 



In order to solve these equations we search for a solu- 
tion in form of an expansion of regular Bessel functions, 
so that we write 



*(f) = *(r,0)= ^ aiJi{kr)e'' 



(3) 



l = -L 



This expansion is obviously a solution of the first equa- 
tion, so that we have to solve it under the boundary con- 
ditions Eq. (y). On the billiard boundary, which may be 
parameterized by the curvilinear abscissa s, the expan- 
sion Eq. (H) becomes 



$(s) = *(r(s),0(s)) = ^ aiJi{kr{s))e' 



■W(s) 



(4) 



l=-L 



which is a periodic function of s. This periodic function 
may be decomposed in a Fourier series whose coefficients 
are given by 



Cnik) 



2n 



ds^{s)e' 



-2in7rs/ C 



J2 ^i^nAk) (5) 



CnAk)^ dse-^'^^'/^Mkr{s))e'''^'^ (6) 



with 



To satisfy the boundary condition Eq. (0), one may 
impose the equivalent conditions that all Fourier coeffi- 
cients Cn{k) are zero, at least for —L < n < L, so that 
we are left with the following linear system in the af. 



J2 aiCn,iik)=0,-L<n<L 



(7) 



l=-L 



For this homogeneous system to have a non trivial so- 
lution, one has the following set of equations: 



D{k) = det[CnAk)] = Q 



(8) 



Thus one is left with the problem of constructing the 
matrix Cn,i{k) E^iid finding the zeroes of its determinant 
D(k). Before going further, we would like to point out 
the advantages of this method as compared to the well 
known collocation method. In our method, one may vary 
independently the number of boundary points used to 
evaluate the integrals Cn,i{k) and the number of partial 
waves used in the expansion, whereas this is not so in 
usual boundary methods. (For other methods we refer 
the reader e.g. to Refs. pjl^). Furthermore, it appears 
clearly that if the billiard is close to a circle, the matrix 
is nearly diagonal so that its determinant is easy to com- 
pute numerically, whereas this is not true for plane wave 
decompositions. Finally, in contrast with the Green's 
method, one may always deal with real matrices by de- 
composing on sine and cosine rather than on exponen- 
tials (for time reversal invariant systems). However, one 
should note that this method does not work for billiards 
whose boundary consists of several distinct curves (for 
instance Sinai' billiards) . 



B. Practical 



Therewith one tabulates the function D{k) in an inter- 
val of k such that the number of partial waves needed is 
constant in this interval: 
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(10) 



(11) 



In this equation, Rmax is the greatest distance of the 
boundary from the origin, in our case Rmax = R + a. 
The number of coefficients a„ may then be taken as 



N = L/2 + dN 



fl2) 



where dN typically ranges from to 3 or 4. Usually, 
the positions of the eigenvalues depend very little on this 
parameter. 
The integrals 

Cn,m{k) = / ds sin{mTs/C)J2m{kr{s))sin{2m6{s)) 
Jo 

(13) 

are evaluated using points regularly spaced along the 
boundary. As soon as their spacing is such that the 
fastest varying phase in Eq. ([l3| ) changes by less than 
TT in a given step, the evaluation of these integrals is ac- 
curate enough to give the position of the zeroes of D{k). 
In other words, the step As used in evaluating Eq. ( p^ ) 
should be such that the following condition is verified: 
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(14) 



where dmin is the smallest distance of the boundary from 
the origin, in our case R. With these ingredients, the pre- 
cise computation of the levels of a rather regular billiards 
is fast, and does not require powerful computers. Typi- 
cally, for the stadium, the computation of the 1000 first 
levels takes a few hours run on a personal computer, the 
dimension of the matrices used being smaller than 100. 
The precision obtained is much better than 1/lOOth of 
the average spacing, which is enough for our purpose. 



In practice, the application of the above algorithm de- 
pends on the problem one has to solve. For instance, 
a particular choice of the origin of the coordinates may 
simplify notably the evaluation of the matrix: For the 
stadium, a good choice would be the symmetry center 
of the billiard, so that one may separate easily different 
symmetry classes. With this choice, the expansion given 
in Eq. (pi) become for odd-odd symmetry: 



N 



*(f) = *(r,6l) == ^a„J2„(fcr)sm(2n6l) 



(9) 



IV. SPECTRAL FLUCTUATIONS 

In the following section we want to discuss the sta- 
tistical properties of the data obtained for the 7 = 1 
stadium billiard as described in Sees. H and III. After a 
short summary of the general concepts of analyzing spec- 
tral fluctuations in Sec. [V A , we present the results for 
the comp arison of the meas ured and calculated data in 



Sec. [VB. Finally, in Sec. IV C we compare the influ 



n=l 



ence of a special family of orbits in different Bunimovich 
stadium billiards. 



A. Theoretical background 

From the measured resp. numerically simulated 
eigenvalue sequences ("stick spectrum") the spectral 
level density p{k) —'Y^^5{k — ki) is calculated (fc is the 
wave number, fc — 27r/c • /) and a staircase function 
N{k) = J p{k')dk' is constructed which fluctuates around 
a smoothly varying part, defined as the average of A^(fc). 
Usually this smooth part 7V*™°''*''(fc) is related to the vol- 
ume of the classical energy-allowed phase-space. For the 
2D billiard at hand with Dirichlet boundary conditions 
it is given by the Weyl-formula |l8|Jl9|] 
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(15) 



where A is the area of the billiard and C its perime- 
ter. The constant term takes curvature and cor- 
ner contributions into account. Higher order terms 
pO| are not relevant for the present analysis. The 
remaining fluctuating part of the staircase function 
Nf^^<={k) =N{k)-N^''y^{k) oscillates around zero. 
While Eq. (ITq) does not contain any information regard- 
ing the character of the underlying classical dynamics of 
the system, the fluctuating part does. 

In order to perform a statistical analysis of the given 
eigenvalue sequence independently from the special size 
of the billiard, the measured resp. calculated spectrum is 
first unfolded pi|, i.e. the average spacing between ad- 
jacent eigenmodes is normalized to one, using Eq. (15). 
This proper normalization of the spacings of the eigen- 
modes then leads to the nearest neighbour spacing distri- 
bution P(s), from now on called NND, the probability of 
a certain spacing s between two adjacent unfolded eigen- 
frequencies. To avoid effects arising from the bining of 
the distribution, we employ here the cumulative spacing 
distribution I{s) = J P{s)ds. 

To uncover correlations between nonadjacent resonan- 
ces, one has to use a statistical test which is sensitive on 
larger scales. As an example we use the number variance 
S^ originally introduced by Dyson and Mehta for studies 
of equivalent fluctuations of nuclear spectra |§,^ . The 
S^-statistics describes the average variance of a number 
of levels n{L) in a given interval of length L, measured 
in terms of the mean level spacing, around the mean for 
this interval, which is due to the unfolding equal to L, 

Y.\L) = l^{n{L)-{n{L))Lf^ ={n^{L))L-L\ (16) 

Furthermore, to characterize the degree of chaoticity 
in the system, the spectra are analyzed in terms of a sta- 
tistical description introduced in a model of Berry and 
Robnik p3| which interpolates between the two limit- 
ing cases of pure Poissonain and pure GOE behavior for 
a classical regular or chaotic system, respectively. The 
model introduces a mixing parameter q which is directly 
related to the relative chaotic fraction of the invariant Li- 
ouville measure of the underlying classical phase space in 



which the motion takes place (5 = stands for a regular 
and g = 1 for a chaotic system). 

The chaotic features of a classical billiard system are 
characterized by the behavior of the orbits of the prop- 
agating point-like particle. The quantum mechanical 
analogue does not know orbits anymore but only eigen- 
states, i.e. wave functions and corresponding eigenener- 
gies. Thus, the individual and collective features of the 
eigenstates must reflect the behavior of the classical or- 
bits. The semiclassical theory of Gutzwiller m assumes 
that a chaotic system is fully determined through the 
complete set of its periodic orbits. The influence of the 
isolated periodic orbits of a billiard is most instructively 
displayed in the Fourier transformed (FT) spectrum of 
the eigenvalue density p-'''"^(fc) = dN^^'^'^{k)/dk, i.e. 



p^'"=(x) 



[pik)-p'^^y\k)]dk, (17) 



with [fc. 

which the data are taken 



mm,kmax] being the wave number interval in 



B. Results 



In this section we present the results of the analysis 



of our data using the techniques described in Sec. IV A 



We start with a direct one-to-one comparison of both, 
experimentally and numerically obtained data sets, and 
perform first a comparison with the sequences of the 
unfolded eigenvalues. To obtain the unfolded eigen- 
values, we fit the Weyl-formula, Eq. ([l5|), onto our 
measured resp. numerically simulated spectral stair- 
case N{k). Doing this, one obtains the following pa- 
rameter of Eq. ( |l5| ) for the measured data: an area 
Afit.exp — (710.52±4.12) cm^ and a perimeter Cfn^exp = 
(113.68 ± 2.16) cm, which are very close to the design 
values Adesign = 714.15 cm^ and C design = 111.41 cm. 
Within the given uncertainties the respective coefficients 
A and C agree fairly well. From the numerical data which 
were obtained from a 7 = 1 stadium billiard with area 
A — An the coefficient C can be calculated. The fit of 
the Weyl-formula to the numerical data yields values for 
coefficients A and C within a few per mille of the calcu- 
lated coefficients. From this we conclude that due to the 
not well known mechanical imperfections of the billiard 
the uncertainties in the average properties of the spec- 
trum expressed through the coefficients A and C of the 
Weyl-formula are larger in the experiment than in the 
numerical simulations. 

With the unfolding indicated above the different ge- 
ometrical properties of the real and ideal billiards are 
removed and a direct comparison of the first 770 eigen- 
values becomes possible. Computing the difference of the 
experimental and the numerical data, e. 
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plotting this difference over the unfolded experimental 
eigenvalues, the upper part of Fig. 2 is obtained. The 



curve fluctuates around zero, which is an indication for 
complete spectra, i.e. no mode is missing. The appearing 
oscillations in the curve are random and reflect the fact 
of how accurate a certain eigenfrequency could be deter- 
mined experimentally. If only one eigenvalue in the mea- 
sured sequence is artificially removed, the curve shows a 
clearly observable step of height 1 at the position where 
the eigenvalue was dropped. This is demonstrated in the 
lower part of Fig. 2, where the 427*'' mode has been re- 
moved in the experimental spectrum. More information 
can not be extracted from this direct comparison since 
systematic errors and imperfection have been removed 
through the unfolding procedure. Therefore we want to 
concentrate in the following on the statistical properties 
of the presented complete sequences of data. 
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FIG. 2. In the upper part the difference of the unfolded 



eigenvalues 



is plotted over the unfolded eigen- 



frequencies eexp of the experimental data. In the lower part 



the same curve is shown, but with the 427 
removed from the experimental spectrum. 



mode artificially 



In the next step, we extract the fluctuating part of the 
staircase function iV/'"'=(fc) = N{k) - N^^v^ (fc) from the 
measured resp. calculated eigenvalue sequences (Fig. 3). 
As expected oscillations around zero are seen, also indi- 
cating no missing modes. 

In the upper part of Fig. 3 a strong enhancement of 
the amplitude of the fluctuations as a function of wave 
number fc - a characteristic feature of regular systems - is 
observed as well as a periodic gross structure with spac- 
ing of 15.7 m-i (750 MHz). This behavior of 7V^'"=(fc) 
is caused by the well known family of marginally stable 
periodic orbits, which bounce between the two straight 
segments of the billiard, the so-called bouncing ball or- 
bits (bbo), with length 2r. The cumulative level density 
N{k) shows therefore periodic oscillation with a flxed pe- 
riod of 27r/2i? = 15.7 m~^ around the value given by the 
Weyl-formula. These observations can be described |2J] 
by the semiclassical expression of the contribution of the 
bbo to the spectrum which reads 
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(18) 



with X = (kr)/Tr. This formula (for a generalization of it 
in three dimensions see [[L2[) reproduces the mean behav- 
ior of the experimental data as shown in the upper part 
of Fig. 3. After subtraction of this smooth correction in 
addition to the expression of Weyl the proper fluctuat- 
ing part of the level density is obtained, which is plotted 
in the lower part of Fig. 3. Naturally the same result is 
obtained for the data from the numerical calculations. 
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FIG. 3. The histogram in the upper part displays the fluc- 
tuating part N{k) — N^''^\k) of the staircase function. The 
full line shows the semiclassical prediction for the bouncing 
ball orbits (bbo) according to Eq. (|l8[). The lower part shows 
N^'^^ik) - N''''°{k), the fluctuating part after subtracting the 
contribution of the bbo. 



To determine the degree of chaoticity of the investi- 
gated 7 = 1 stadium billiard we next calculate the cu- 
mulative nearest neighbour spacing distribution I{s) for 
the fluctuating part of the staircase function corrected 
by the Weyl and the bbo terms. In Fig. 4 the dashed 
curve shows the density for a certain spacing s of two 
adjacent unfolded eigenmodes for the experimental and 
the numerical data. Beside the data also the two lim- 
iting cases, the Poisson- and the GOE-distribution, are 
displayed. Using the ansatz of Berry-Robnik one obtains 
a mixing parameter q = 0.97 ±0.01 for the experimental 
and q = 0.98 ± 0.02 for the numerical data. As stated 
in PI the Bunimovich stadium billiard should be fully 
chaotic, which is expressed through q being very close to 
unity within the uncertainties of the fltting procedure. 
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FIG. 4. Cumulative nearest neighbour spacing distribution 
for the 7 = 1 stadium billiards. The dashed curve corre- 
sponds to the experimental (left side) respectively the numer- 
ical (right side) data. Also the two limiting cases (Poisson and 
GOE) are displayed. The contributions of the bbo are already 
extracted, so that the data show the predicted GOE-behavior. 

For the E^-statistics which measures long-range corre- 
lations the effects of the bbo's are strikingly visible, see 
Fig. 5. Their presence influences the rigidity of the spec- 
trum for large values of length L. A proper handling of 
these orbits as described for the cumulative NND changes 
the form of the curve towards the expected GOE-like be- 
havior. In the upper part of Fig. 5 the experimental 
data (□) and the numerical data (A) are displayed for 
the E^-statistics before removing the contribution of the 
bbo's. Comparing the E^-statistics of the numerical and 
the experimental data, no significant difference between 
both data sets can be seen. In contrast to this a small 
difference between both is extracted when the bbo con- 
tribution is removed (lower part of Fig. 5). However, in 
both cases the E^-statistics follow very closely the GOE 
prediction up to L w 3.5, where the distribution satu- 
rates. This is in a very good agreement with the theo- 
retical predicted value according to |2^ of L„iax ~ 3.8. 
This length L — L„iax also refers to the shortest periodic 
orbit, which is in the present billiard the bbo. Why is 
E^(L) different for L > L,nax for the two data sets? This 
might be due to the fact that the fabricated microwave 
cavity is a real system with small but existing mechanical 
imperfections, e.g. a slight non-parallelity of the straight 
segments of the cavity, which should have its effect on 
the periodic orbits. 

After comparing the statistical measures of the two 
investigated 7 = 1 stadium billiards, we now consider 
their periodic orbits which are given through the Fourier 
transformed spectrum of the eigenvalue density ^•'''"'^(fc), 
see Eq. (|^. The lengths of the classical periodic or- 
bits (po) correspond to the positions and their stability 
roughly spoken to the height of the peaks in the spec- 
trum. In Fig. 6 the mod-squared of the Fourier trans- 
formed oi p^^'^'^{k) for the experimental (upper part) and 
the numerical (lower part) data are compared. The nu- 
merical length spectrum is scaled to the experimental one 
in such a way that the lowest bbo occurs at the same loca- 
tion in the two Fourier spectra. Then the bbo's, because 
of their dominating nature in the spectra, were removed 



p6[ . Now a direct comparison between experiment and 
simulation becomes possible. It can be seen in Fig. 6 that 
beside the positions of the peaks (lengths of the po) also 
their heights (stability of the po) are almost identical in 
both spectra. 




FIG. 5. E^-statistics of the experimental data set (□) com- 
pared with the numerical data (A). In the upper part the 
contribution of the bbo is not extracted, in the lower part it 
is. The error bars of the data points are a measure of the sta- 
tistical fluctations within the given ensemble of eigenvalues 
in an interval of length L (see Eq. (|l6|)). Note, the predicted 
saturation observed in both curves happens at L = Lmax ~ 4. 




FIG. 6. Fourier transform of the fluctuating part of the 
level density p(fc) -p^"=^'(fc) -p''''°(A;). The experimental and 
the numerical results are displayed as mirror images. Rem- 
nants of the bbo's at a:: = 0.4 m, at 0.8 m and to a lesser 
extent at 1.2 m are still visible. 



Let us finally in this subsection return to the bounc- 
ing ball orbits. A detailed analysis of those in the 
experimental billiard, which have lengths x — n ■ R 
[n = 2, 4, 6, . . .) shows that the distance R between the 
straight lines of the microwave resonator is not exactly 
R — 2Q cm as specified in the construction but rather 
R — (19.92 ± 0.05) cm (average over the first six recur- 



rences of the bbo's). That means that the two straight 
segments of the bihiard are not exactly parallel. The 
reasons for this deviation {AR/R = 4 x 10^^) from the 
design value are on the one hand a finite mechanical tol- 
erance during the fabrication (« 3 x 10^'^) and on the 
other hand effects coming from thermal contractions ||2^] 
at low temperatures (« 1 x 10"'^). This small mechani- 
cal imperfection is also the reason for the higher satura- 
tion level of the S^-statistics in the experimental data in 
comparison to the numerical one after removing the bbo 
contribution (lower part of Fig. 5). 



C. Influence of the bouncing ball orbits in different 
Bunimovich stadiums 

In this section we discuss finally the influence of the 
bouncing ball orbits in the length spectra and in the 
level statistics of different Bunimovich stadium billiards. 
Therefore we compare the 7=1 stadium discussed so far 
with an also measured superconducting 7 = 1.8 stadium 
billiard (a = 36 cm and R = 20 cm), already presented 
in g. For this comparison we use the first 1060 eigenfre- 
quencies up to a frequency of 18 GHz. 
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FIG. 7. Experimental length spectra for the 7 = 1 (upper 
part) and the 7 = 1.8 (lower part) stadium billiard. In both 
cases the bbo are not extracted. Obviously the spectra reveal 
different amplitudes at the positions of the bbo, which are 
marked by the dashed lines {x = n ■ R, n = 2, 4, 6 . . .). 

The interesting features for comparing the two stadi- 
ums are the amplitudes of the peaks for the bbo and not 
the peaks of the unstable po's, since they naturally have 
to differ. In Fig. 7 the length spectra of the 7 = 1 and 
the 7 = 1.8 stadium billiards are shown. It is clearly 
visible that the peaks of the bbo's have different heights 
or amplitudes in the two investigated stadiums. Because 
the height of the peaks corresponds to the stability of the 
periodic orbits, these differences can be easily explained. 
For the 7 = 1.8 stadium the geometrical area on which 
the bbo's exist is almost twice the area of the 7 = 1 



billiard. This correlation between the area of the rect- 
angular part of the billiards with radius R fixed and the 
amplitude of the peaks for the bbo in the FT has been 
first pointed out in (28|. This effect is also expressed 
through the ratio a/R in Eq. (18). 




FIG. 8. Comparison of the E^-statistics of the 7 = 1.8 sta- 
dium (o) before (upper curve) and after (lower curve) the 
extraction of the bbo's contribution, as well as for the 7 = 1 
stadium (□). 

The effects of the bbo's are furthermore visible in 
the statistical measures, as pointed out in the previ- 
ous section. For the nearest neighbour spacing distri- 
bution we have found a Berry-Robnik mixing parame- 
ter q — 0.97 ± 0.02 for the experimental 7=1 billiard 
before and q = 0.98 ± 0.01 after removing the bbo con- 
tribution. But for the j — 1.8 stadium, presented in 
Q, the situation differs. There we have a mixing pa- 
rameter q = 0.87 ± 0.03 for the billiard with bbo's and 
q ~ 0.97 ± 0.02 without them. Using the S^-statistics 
to investigate the long range correlations, the effect of 
the bbo's for stadiums of different 7 values becomes even 
more obvious. The S'^-statistics for the 7 = 1.8 stadium, 
including bbo's, increases more strongly, see upper part 
of Fig. 8 (o), than for the 7=1 billiard in the lower part 
of the same Figure (□). After removing the bbo contribu- 
tion one notices saturation sets in at the same L„iax ~ 4 
for both stadiums. In other words, the shortest po is 
the same in both stadiums as it should be. The results 
which have been found are thus in excellent agreement 
as predicted in |^,^ which state that the influence of 
the bbo's is in the 7=1 stadium billiard is very small, 
so that it can be characterized as being the most chaotic 
stadium. 



V. CONCLUSION 

In this paper we have shown two different methods, 
an experimental and a numerical one, to obtain eigen- 
value sequences of a 7 = 1 Bunimovich stadium billiard. 
A direct comparison of these two data sets clearly re- 
veals that both sets are complete and that no eigenmode 



is missing. Informations, e.g. accuracy of the measured 
resp. simulated eigenvalues can not be obtained from this 
test. Therefore we analyzed the statistical properties of 
the billiards. Using the cumulative nearest neighbour 
spacing distribution I{s) both data sets show the same 
predicted GOE behavior after removing the contribution 
of the bouncing ball orbits. The same result is obtained 
from the S~-statistics, by which we investigated the long 
range correlations of the eigenvalues. Only the saturation 
level of the experimental data at large intervals is some- 
what above of the level of the numerical data, because 
the bbo's can not be removed exactly. The reason can be 
found in the length spectra of the data which uncovers 
that the lengths of the bbo's in the experiment slightly 
deviate from their expected values due to the design of 
the cavity. The cavity contracts at low temperatures and 
its mechanical fabrication has only been possible within 
certain tolerances. On the other hand we have shown 
that such small imperfections of the real system do not 
have any influence on the results of the statistical anal- 
ysis. Furthermore, we have investigated the influence of 
the bbo in two different stadium billiards. They have a 
clearly different effect on the respective statistical mea- 
sures, confirming in particular that the 7=1 stadium is 
even more chaotic than the 7 = 1.8 one. 

The comparison of the experimental and the numeri- 
cal eigenvalues was performed to show how accurate the 
results of the statistical analysis depend on the given 
method of providing eigenvalues. To obtain informations 
about simple two-dimensional billiards, like the presented 
7 = 1 stadium billiard, numerical calculations some- 
times have an advantage, e.g. if one is interested only 
in eigenvalues. To simulate experiments involving eigen- 
functions like in [|l3|,0| is a different matter. Further- 
more, for problems where one is interested in billiards 
with scatters inside, billiards with fractal boundaries or 
three-dimensional billiards, etc. the experiment clearly 
offers a very convenient way to obtain large sets of eigen- 
values quickly. 
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